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ABSTRACT 



We report the observations of PG 1553+113 during the first ~200days of 
Fermi Gamma-ray Space Telescope science operations, from 4 August 2008 to 
22 February 2009 (MJD 54682.7-54884.2). This is the first detailed study of 
PG 1553+113 in the GeV gamma-ray regime and it allows us to fill a gap of three 
decades in energy in its spectral energy distribution. We find PG 1553+113 to 
be a steady source with a hard spectrum that is best fit by a simple power-law in 
the Fermi energy band. We combine the Fermi data with archival radio, optical, 
X-ray and very high energy (VHE) gamma-ray data to model its broadband 
spectral energy distribution and find that a simple, one-zone synchrotron self- 
Compton model provides a reasonable fit. PG 1553+113 has the softest VHE 
spectrum of all sources detected in that regime and, out of those with significant 
detections across the Fermi energy bandpass so far, the hardest spectrum in that 
energy regime. Thus, it has the largest spectral break of any gamma-ray source 
studied to date, which could be due to the absorption of the intrinsic gamma-ray 
spectrum by the extragalactic background light (EBL). Assuming this to be the 
case, we selected a model with a low level of EBL and used it to absorb the power- 
law spectrum from PG 1553+113 measured with Fermi (200 MeV- 157 GeV) to 
find the redshift which gave the best fit to the measured VHE data (90 GeV - 
1.1 TeV) for this parameterisation of the EBL. We show that this redshift can be 
considered an upper limit on the distance to PG 1553+113. 

Subject headings: Gamma rays: observations — BL Lacertae objects: individual 
(PG 1553+113) 



Introduction 



PG 1553+113 is a high-frequency peaked BL Lacertae object (HBL; iFalomo et al.lll994 



Beckmann et al.ll2002l ) whose redshift remains unknown despite continued efforts. Like all 
BL Lacs, we find its spectral energy distribution (SED) to have a double-peaked shape 



51 Consorzio Interuniversitario per la Fisica Spaziale (CIFS), 1-10133 Torino, Italy 
52 Dipartimento di Fisica, Universita di Roma "Tor Vergata", F00133 Roma, Italy 
53 School of Pure and Applied Natural Sciences, University of Kalmar, SE-391 82 Kalmar, Sweden 
* Corresponding authors: D. Horan, deirdre@in2p3.fr; D. Sanchez, sanchez@poly.in2p3.fr 
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(in vFv representation) that can thus be parameterised with four characteristic slopes. It 
has been detected from radio through hard X-r ays and also in the very high energy (VHE ; 
E > 100 GeV) regime up to energies above 1 TeV ( Aharonian et al]|2006b ; Albert et al. 2007 ). 
With these data three of the four components of its SED were sampled, namely, the rising 
and falling edges of the low-energy peak (~10~ 6 eV - 30 keV) and the falling edge of the high- 
energy peak (~90 GeV - 1 TeV). We report here the first detailed analysis of the rising high- 
energy portion and, crucially, the high-energy peak, of the PG 1553+113 SED. These data, 
from observa tions made by the Fermi Gamma-ray Space Telescope Large Area Telescope 
(Fermi LAT; lAtwood et al.ll2009l ) during its first ~200 days of operation, are combined with 
data from other wavebands to construct and model the broadband SED of PG 1553+113 in 
Section SU 



Discovered and classified as a BL Lacertae object (BL Lac) by iGreen et al.l (119861 ). 
PG 1553+113 is located at R.A. of a J2000 = 15h55m43.04s and declination of 5j2ooo = 
+lldllm24.4s in the constellation of Serpens Caput. The logarithmic ratio of its 5 GHz radio 
flux, F 5 ghz, to its 2keV X-ray flux, Fikev x has been found to range from \og(F 2 keV / F^ghz) 
= -4.99 to -3.88 (lOsterman et alihood : Sector et aDl2003l ). The high value of this ratio 
places PG 1553+113 at times among the most extre me of the HBLs; a BL Lac is classi- 
fied as extreme when it has \og(F 2 k e v / F^ghz) > -4.5 (IRector et al.ll2003l ). A number of the 
TeV blazars have, at one time or another, exhibited fluxes that place them in this extreme 
category (e.g., 1ES 0229+200, H 1426+428, 1ES 1959+650: IRector et allhoO^ . 



In the radio band, PG 1553+113 has been detected at different mean flux levels. Its 



4.8 GHz flux, for example, has ranged from 180 to 675 mJy ( 


Bennett et al.lll986: 


Gregory & Condon 


1991: 


Becker et al. 


1991 




Perlman et al. 


2005: 


Osterman et al. 


2006). Its flux between 4.8 



and 14.5 GHz was fo und to be variable on mon th timescales during the observations of 
Perlman et al.l (120051 ) and lOsterman et al.l (120061 ). VLB A observations have resolved a jet 
extending at least 20 pc to the northeast of PG 1553+113 (IRector et al.ll2003l ). No evidence 
for superluminal motion has been reported in the literature to date; multi-epoch VLBA 
monitoring has commenced recentljQ. 



PG 1553+113 is a brigh t optical source with ^/-band magnitude of V Q ~ 14 (IFalomo fc Treves 



19901 : lOsterman et al.ll2006l ). Observations taken between 1986 and 1991 with the ESO tele- 
scopes found its spectral index, q fl1 to remai n almost constant (a ~ — 1) and its magnitude 
to vary by AV n _ = 1.4 (IFalomo et al.lll994j ). Low levels of optical variability were seen by 



Osterman et al.l (120061 ) during their 2003 observing campaign. Limits on the magnitude of 
its host galaxy will be discussed later. 



1 http://web. whittier.edu/gpiner/research/index. htm 
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PG 1553+113 is also a bright X-ray source that has been observed by most X-ray mis- 
sions (Einstein, ROSAT, ASCA, BeppoSAX, RXTE, XMM-Newton, Swift and Suzaku). Al- 
though it has been detected at a number of different flux levels by these observatories, no 
evidence for strong or fast (sub-hour) flux variability has been observe d at X-ray energies 
(IReimer et al.ll2008l ). The Suzaku observations performed in July 2006 (IReimer et al.ll2008l ) 
provide the highest energy X-ray measurement so far obtained for PG 1553+113 at ~ 30keV 
with a 10-30 keV flux of 1.35 x 10 -11 erg cm -2 s _1 . No evidence for spectral hardening up to 
these energies was found in these data, indicating that all of the X-ray emission detected was 
due to synchrotron emission. The 2-10 keV fluxes measured by different X-ray missions are 
shown in Table [R the highest being 6.9 x 10~ n erg cm -2 s -1 in October 2006 (Swift-XRT; 



Tramacere et al. The X-ray flux w as found to doub l e over a period of 10 days during 



the 3- week RXTE observing campaign of lOsterman et al.l (120061 ). Despite the variations of 
a factor of approximately 20 in the 2-10 keV X-ray flux, the measured spectral properties of 
PG 1553+113 at these energies, also listed in Table [TJ have not changed significantly over the 
course of the X-ray observations. It exhibits spectral curvature that can be well- described 
with either a broken power law or a log-parabolic shape. In the Suzaku data this curva- 
ture extends into the hard X-ray band (< 30keV) and steepening of the spectrum above 
~ 10 keV beyond that predic ted by either the broken power-law or log-parabolic model is 
observed (IReimer et al.l 120081 ) . 

In the high energy (HE; E ~ 100 MeV - 100 GeV) gamma-ray regime, PG 1553+113 
was not detected by EGRET. An upper limit of F EGRET < 9.97 x 10~ 8 cm" 2 s" 1 abov e 



100 M eV was derived based on the summed exposure of cycles 1,2,3 and 4 (IHartman et al. 



19991 ). At higher energies, in the VHE regime, PG 1553+ 113 is a confirmed gamma-ray 



emitter. First detected by H .E.S.S. (jAharonian et al.ll2006bl ) and subsequently confirmed by 
MAGIC (Ulbert et alil2007h . PG 1553+113 has a flux that is approximately 3% that of the 
Crab Nebula at these energies. The combined 2005-2006 H.E.S.S. data fo llow a power law 
with photon index of Py#£=4.46 ± 0.34 between ~225 GeV and 1.3 TeV (jAharonian et al. 



20081 ). The MAGIC 2005-2006 spectra are consistent with this having a power-law photon 
index of T VHE =A.21 ± 0.25 between ~90 and 500 GeV (Ulbert et al.ll2007h . No evidence for 
variability of the spectral index, to within the measurement uncertainties, is seen in the 
VHE measurements. The H.E.S.S. and MAGIC integral flux levels are consistent from April 
to August 2005 and from April to July 2006 with a mean flux of Ivhe (E> 200 GeV) = 
2.0 + 0.8 x 10~ n cm~ 2 s -1 . Although the spectrum remained unchanged, the flux detected 
by MAGIC between January and April 2006 was lower than the preceding and proceeding 
measurements at Ivhe (E > 200 GeV) = 0.6 + 0.2 x 10~ n cm -2 s -1 . When the systematic 
uncertainties on the fluxes are taken into account, this flux is marginally inconsistent with 
the VHE fluxes detected during the H.E.S.S. and other MAGIC observations suggesting 
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that the PG 1553+113 flux varied by up to a factor of three on monthly timescales in 2006. 
PG 1553+113 has the steepest spectrum of all of the sources detected in the VHE regime, 
which makes it a promising target for the Fermi LAT because extrapolating down to the 
Fermi energy range would make this an extremely bright source unless a dramatic spectral 
break occurs at energies below ~100 GeV. EGRET's non-detection could be interpreted in 
this context or as the result of PG 1553+113 being in a lower emission state at that time than 
during the VHE observat ions. PG 1553+113 is in the Fermi LAT bright AGN source list 



(LBAS; lAbdo et al.ll2009al ) . the high-confidence AGN associations from the first three months 
of Fermi data. Listed as OFGL J1555.8+1110 with a flux of I LB as (E > 100 MeV) = 8.0 ± 1.0 
x 10~ 8 cm' 2 s' 1 and a photon index of T^bas — 1-70 ± 0.06, it has one of the hardest spectra 
of the 106 AGN that comprise this list. Indeed, if only the AGN with a significant detection 
across the entire Fermi bandpass are considered, its spectrum is the hardest of those in the 
LB AS. This combination of a very soft VHE spectrum and a very hard HE spectrum means 
that PG 1553+113 has a significant spectral break in the gamma-ray regime. 

Despite significant efforts, the redshift of PG 1553+113 remains unknown. Its measure- 
ment is of great interest both for a better understanding of its SED, in particular since it 
is an extreme BL Lac and therefore has a very hard synchrotron emission spectrum, and 
also for studying EBL effects; some redshift estimat es make it the most distant VHE source 



detected to date at z <0.78 (ISbarufatti et al.ll2005l ). which, if shown to be the case, could 



imply an EBL density close to the minimum allowed by galaxy counts. The results of red- 
shift studies undertaken to date are summarized in Section [3] and the archival gamma-ray 
data together with the Fermi data from PG 1553+113 are used to place constraints on its 
redshift. 



2. Fermi Observations of PG 1553+113 



The characteristics and performance of t he LAT Instrument on the Fermi Gamma- 
ray Space Telescope are described in detail by lAtwood et al.l (120091 ) . Presented here is the 
analysis of the Fermi LAT data from a region of 10° radius centred on PG 1553+113 from 4 
August 2008 to 22 February 2009 (MJD 54682. 7-54884.2 jj These data were analyzed using 
the standard Fermi analysis software^. Events with zenith angles < 105° were selected from 
the so-called "diffuse class" , those events having the highest probability of being a photon. 
Only events with energies greater than 200 MeV were used in this analysis. The extragalactic 



2 Mission Elapsed Time 239,557,414 to 256,966,310 seconds. 

3 ScienceTools v9rl0 with instrument response function (IRF) P6_V3_DIFFUSE ( Atwood et al.ll2009h . 
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diffuse gamma-ray emission together with the residual instrumental background was modeled 
as a simple power-law while the galactic diffuse was modeled with GALPRoM 

([Strong et al.l 



2004b] , IStrong et al.ll2004a[ ) 



There are two other Fermi sources within 10° of PG 1553+113. These sources, lying 
at angular separations of 1.8° and 5.5° from PG 1553+113 were modeled in our analysis so 
that they could be subtracted out along with the backgrounds described above. The nearest 
source to PG 1553+113, OFGL J1553. 4+1255, is located at R.A. of a J20 oo = 15h53m28.2s 
and declination of <5/2ooo — +12d55m20.3s and is thus spatially coincident with the quasar 
QSO B1551+1305 (z=1.29). It is detected up to approximately 10 GeV with an integral 
flux (E >200MeV) of 5.67 + 0.38 x lO^cm^s" 1 and photon index of 2.26 + 0.05. There 
is evidence for variability in its photon index and flux but neither are correlated with the 
signal detected from PG 1553+113. The LAT has a point spread function (PSF), # 68 , that 
decreases with increasing energjIE For energies above 360 MeV, the PSF of the LAT is 
smaller than 1.8°, the angular separation between PG 1553+113 and OFGL J1553. 4+1255. 
Therefore, for almost all of the energy range over which the PG 1553+113 data are analyzed 
here (E = 200 MeV - 400 GeV) , the PSF of the LAT is such that the PG 1553+113 data are 
not significantly contaminated by the signal from OFGL J1553. 4+1255. To ensure that this 
was the case, we performed an analysis of the PG 1553+113 data with energy E > 400 MeV. 
The results obtained from this analysis are entirely consistent with those obtained when the 
full energy range (E > 200 MeV) is considered. 

The other source within 10° of PG 1553+113 is not in the LBAS. It is located at a R.A. 

of 

«J2ooo = 16h07m40.5s and declination of <5j2ooo — +15d50m36.2s and was detected up to 
approximately 10 GeV with an integral flux (E >200MeV) of 2.06 + 0.28 x lO^cm^s -1 
and photon index of 2.31 ±0.11. Its angular separation of 5.5° from PG 1553+113 is well in 
excess of the width of the Fermi PSF over the full energy range. 

In the Fermi data analysed here, PG 1553+113 was detected with a statistical signifi- 
cance of approximately 49a (gtlike calculated a test significanc^l of 2407) and an integral 
flux /(£>200MeV) of 5.00 + 0.31 x 10" 8 cm~ 2 s _1 . The most energetic photon in the 
Fermi data is at 157 GeV (at an angular separation of 0.04° from the source location, i.e., 



4 The background file mapcube_54_59Xvarh7S . f its was used. 

5 The accurate parameterisation of the LAT PSF to be used for science analysis is described by the 
instrument IRF. The following simplified, acceptance-averaged approximation for the 68% containment angle 

might be useful as an illustration of the PSF energy dependence: 9q$ ~ ^/(0.8° x i?Q°y ) 2 + (0.07 ) 2 . 

6 The test significance is defined as - 2 1n(A), where A is the likelihood ratio for the null hypothesis and the 
assumed source model; see lAbdo et al.l (|2009bl ) for a full description. 
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well within #68 for that energy). The Fermi data are well described by a power-law such 
that the differential photon flux, F(E), is given by: 



dE u V^o 

where F is the differential flux at energy, E and F is the photon index. For each energy, 
E, the uncertainty contours, called the butterfly, are defined such that the differential flux 
satisfies: 

where cov(F ,r) is the covariance term, returned by the MINUIT minimization and error 
analysis function called by the Fermi likelihood analysis tool, gtlike, and AF, AF and 
Ar are the statistical uncertainties on the F, F and T, respectively, at energy, E. Equation[2] 
reaches a minimum at the decorrelation energy, E dec , where: 

E dec = £ exp ( C °^° r2 r) ) , F dec = F{E dec ) (3) 

For the PG 1553+113 data analyzed here, we find a differential flux of F dec = 2.60 ± 0.18 
x 10~ 9 cm" 2 s" 1 GeV" 1 with E dec = 2AGeV and a photon index of T = 1.68 ±0.03. The 
differential energy spectrum with the butterfly is shown in Figure [TJ Also plotted are the 
fluxes calculated when the data were analyzed in eight independent energy bins fixing the 
photon index in each bin to that derived from the entire dataset, since the preferred fit to this 
was a power-law. These fluxes are listed in Table [2j To estimate the systematic uncertainties 
on the PG 1553+113 flux and photon index, the PG 1553+113 data were re-analyzed using 



two new sets of instrument response functions (IRFs; lAtwood et al.ll2009l ) that were created 
by propagating both extremes of the uncertainties on IRF P6_V3_DIFFUSE. From this analysis, 
we estimate the systematic uncertainties to be at the level of 2% on the photon index 
(r = 1.68±0.03 stat t°;[! 2 S j /st ) and 3% on the differential flux (F dec = 2M±0A8 stat ±° f 8syst x 
10~ 9 cm" 2 s" 1 GeV- 1 ). 

The PG 1553+113 200-day lightcurves for the integral flux, I 2Day (E >200 MeV), and the 
photon index, T 2 Day, with 2-day binning are shown in Figure [2j This timescale was chosen 
so that a detection with a statistical significance of at least 3a was obtained in most bins. 
These are the longest continuous lightcurves ever derived for this source. The x 2 probability 
of the data being fit by a constant are 0.54 and 0.99, respectively, for the integral fluxes and 
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the photon indices. The normalised excess variance (IVaughan et al.ll2003l ) is cf% xs = —2.7 



x 10~ 2 for the flux lightcurve and <J 2 NXS = — 1.6 x 10 -2 for the photon index lightcurvqll 
Both of these calculations suggest that, to within measurement uncertainties, the flux and 
photon index from PG 1553+113 were constant during these observations. Assuming this 
to be the case, an upper limit on the intrinsic variability that could be present in the flux 
lightcurve was derived as follows. An ensemble of lightcurves, each with a different level 
of intrinsic Gaussian noise added, was simulated. For each level of intrinsic variability, Vi, 
1000 lightcurves, each the same length, L, as the measured lightcurve were generated. Each 
flux point in these lightcurves, (I S im)ij ± (^hD ay )j (1 < j < L), was drawn from a Gaussian 
distribution centred on the mean of the measured lightcurve with width of ^(A/| Day )^. + t> 2 , 

where (Al2Day)j is the measurement uncertainty for the jth measured flux. For each of these 
simulated lightcurves, the excess variance was calculated, providing a distribution of 1000 
simulated excess variance measurements for each level of added intrinsic variability. For 
each of these distributions, the value above which 95% of the excess variances lay was found. 
With these data, upper limits were calculated by constructing a Neyman confidence belt. 
It was found that, at the 95% confidence level, the upper limit on the normalised intrinsic 
variance from the PG 1553+113 flux is v < 7.8 x 10~ 2 . The same procedure was applied 
to the photon index lightcurve and an upper limit of v < 1.4 x 10~ 2 at the 95% confidence 
level was derived for the normalised intrinsic variance on the PG 1553+113 photon index. 



3. Using the Gamma-ray Data to Constrain the Redshift of PG 1553+113 



For many years PG 1553+113 was thought to lie at a redshift of z=0.36 (jMiller &: Green 



19831 ). This redsh ift estimate was based o n a spurious emission line in its spectrum, misidenti- 



fied as Lyman-a ( 



^alomo &: Treves! 



spectral features (IFalomo fc Treves 



99 



(^. Subsequent observ a tions have failed to rey eal any 



199CH : IFalomo et al.lll994j ; ICarangelo et all 120031 ). Since 



its detection in the VHE band, a number of dedicated observing campaigns and reanalyses of 
Hubble Space Telescope (HST) images have been undertaken to dete rmine its redshift but 



despite these efforts, the redshift of PG 1553+113 remains unkn own, (ISbarufatti et al.l 12005 



Sbarufatti et al.ll2006l ; iTreves et al.ll2007l ; lAharonian et al.ll2008l ). Constraints, some of them 
contradictory, have been placed on its redshift using a variety of techniques. These fall into 
two categories — spectral and imaging observations in optical to ultraviolet wavebands and 
VHE gamma-ray spectral observations. 



7 It is possible to arrive at a negative value for <J%xs wnen there is low intrinsic variance and/or the 
measurement uncertainties are overestimated. 
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The host galaxy of PG 1553+113 was neither resolved by lHutchings fc Nefi (119921) using 



the H igh-Resolution camera on the Canada France Hawaii Telescope nor by 



Scarpa et al. 



(120001 ) who observed it as part of a survey of 110 BL Lacs using HST. The h ost galaxies of 69 



out of these were resolved by HST, including almost all objects with z < 0.5. ISbarufatti et al. 



( 120051 ) analyzed these HST data and showed that the host galaxy luminosity is encompassed 
in a relatively narrow range for the BL Lacs resolved with HST, thus concluding that BL Lac 
host galaxies typically have absolute magnitude of M R = —22.8. They used this assumption 
to derive lower limits on the redshifts for those objects lying at unknown distances for which 
no host galaxy could be resolved, thus arriving at a lower limit of z > 0.78 for PG 155 3 +113 , 
the brightest such object in the survey. Using a similar strategy, ICarangelo et al.l (120031 ) 
set a lower limit of z > 0.3 on the redshift of PG 1553+ 1 13 us ing preliminary results from 
observations with the ESO 3.6 m telescope. iTreves et al.l (120071 ) reanalyzed the HST image 
of PG 1553+113 and, again making the assumption that its host galaxy is typical of BL 
Lacs, placed a lower limit of z > 0.25 on the redshift. In comparing their results with 
those from VHE observ ations, they concluded that its redshift is in the range z = 0.3 — 0.4. 
Sbarufatti et al.l (120061 ) used the ESO VLT to measure the optical spectra of a number of 
BL Lacs including PG 1553+113. Combining the fact that no spectral features were resolved 
with their knowledge of the sensitivity of the VLT to such features, they were able to derive 
a lower limit of z > 0.09 on the redshift PG 1553+113, under the assumption that its host 
galaxy is typical. 



All of the results above, however, rely on the result of ISbarufatti et al.l (120051 ) . which 
states that BL Lac host galaxies have a very sm all dispersion in their absolute magnitudes. 
The results of a study by lO'Dowd fc Urryl (120051 ) casts doubt on the validity of this assump- 
tion for objects with z > 0.5 because they find evidence for strong evolution in the host 
galaxies of BL Lacs in the redshift range of z — 0.5 — 2.5. 

Gamma-ray data can be used indirectly to estimate redshifts provided that a few 
necessary assumptions are valid. The VHE spectra of extragalactic sources whose red- 
shifts are known can be used to probe the density of the EBL because the high-energy 
gamma rays produce e + e" pairs with the EBL photons thus introdu cing a redshift-dependent 



absorption feature on the spectra observed i n the VHE regime (|Costamante et al.l 120041 ; 
Dwek fc Krennrich 2005 ; Stecker et al. 2007 ; Krennrich et al. 20081). Conversely, in cases 



such as that of PG 1553+113, a firmly established VHE s ource (lAharonian et al.l l2006bl ; 
Albert et al.ll2007l ; lAharonian et al.ll2008l ; lAlbert et al.ll2009l ) lying at an unknown distance, 
the measured spectrum can be used in combination with assumptions about the density 
of the EBL and the i ntrinsic VHE spectru m to put limits on its redshift. By assuming a 



minimal level of EBL (jPrimack et al.ll200ll ) and that the intrinsic spectrum of PG 1553+113 
is unlikely to have a photon index harder than Ti nt = 1.5, the limit from shock accelera- 
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tion models (lAharonian et al.l l2006al ). lAharonian et al.l (l2006bl ) derived an upper limit of 
z < 0.74 on its redshift. We note that none of the Fermi LBAS sources that have significant 
detections across the entire Fermi bandpass have photon indices harder than 1.5, providing 
additional support for the validity of this assumption. When these data were reanalysed and 
combined with subsequ ent H.E.S.S. data, a refined upper lim it of z < 0.6 9 was calculated 
( lAharonian et al.l 120081 ) . A simi lar procedure was adopt ed by lAlbert et al.l (120071 ) resulting 
in an upper limit of z < 0.78. iMazin fc Goebell (120071 ) combined all of the existing VHE 
gamma-ray data from PG 1553+113 and, by assuming that a break in the intrinsic spec- 
trum should be visible if the source lies above z = 0.42, as well as a limit on the hardness 
of the intrinsic spectrum T int = 1.5, placed an upper limit of z < 0.42 on the redshift of 
PG 1553+113. 

The gamma-ray SED measured by Fermi, H.E.S.S. and MAGIC are plotted together 
in Figure [3j The VHE data are non-simultaneo us with those from Fermi and comprise 



two H.E.S.S. datasets and two MAGIC datasets (Aharonian et al. 2006b; Aharonian et al. 



20081 : lAlbert et al.l 120071 : I Albert et al.l 120091 ) . The VHE spectra recorded by H.E.S.S. and 



by MAGIC were consistent during all observations. Although one of the MAGIC datasets 
showed evidence for a change in the mean flux level from PG 1553+113 on month timescales, 
no evidence for strong variability or day-scale (or shorter) VHE flux variability is seen, within 
the measurement uncertainties, in any of the datasets. As described in Section [2] no flux or 
spectral variability was detected in the Fermi dataset. We therefore make the assumption 
here that the source was in a non-flaring state during all of the gamma-ray observations. 
Furthermore, since the highest energy Fermi datapoints, which overlap with the energy range 
covered by MAGIC, show a flux level consistent with that measured by MAGIC, we assume 
that PG 1553+113 was in a similar flux and spectral state at all epochs plotted in Figure [3j 
We note that in both the optical and X-ray regimes, the spectral properties of PG 1553+113 
were not seen to change significantly even when its flux level changed. 

The spectrum measured by Fermi is best-fit by a simple power law. We make the as- 
sumption that any departures from this power-law spectrum up to 1 Te V are dominated by 



abso rption of gamma rays by the EBL and use the parameterisation of iFranceschini et al. 
( 20081 ) . which includes evolutionary effects, to find the level of EBL, and therefore the red- 



shift, which best fits the measured data. By absorbing the extrapolated Fermi spectrum 
with EBL corresponding to redshifts from z = 0.01 — 3.00, in steps of z = 0.01, we find 
that a re dshift of z = 0.75 g ives t he best \ 2 fit to the measured VHE data for the EBL 
model of IFranceschini et al.l (120081 ) . When the Fermi 68% uncertainty-contours, derived 



from Equation [31 are subjected to the same fitting procedure as t he power-law spectrum we 
obtain the statistical error on the redshift determination for the IFranceschini et al.l (120081 ) 
EBL giving us a range from z = 0.70 — 0.79, as illustrated in Figure [3j Since the model 
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of iFranceschini et al.l (120081 ) provides the lowest level for the EBL over the range of interest 
here, the redshift derived from it should be considered an upper limit. Due to the fact that 
we performed this calculation for just one parameterization of the EBL we do not estimate 
the systematic uncertainty associated with absorbing the Fermi spectrum with different EBL 
models and we discuss this further in Section [5j We note that the systematic uncertainty 
on the energy scale, wh i ch is on the order of 10 — 15% for the VHE gamma-ray telescopes 
( lAharonian et al.ll2006d ; lAlbert et al.ll2008l ). was not taken into account in these calculations. 



4. Discussion 



The Fermi data presented here allowed us to derive the longest, continuously sampled 
lightcurves to date for PG 1553+113. Its flux in this energy regime is such that its variability 
can be probed on ~day timescales. The combination of the low duty cycle of VHE instru- 
ments and the weak flux from PG 1553+113 at those energies means that the timescales 
accessible with the LAT for this source are shorter than those accessible at higher energies 
where longer integration times than those available on nightly timescales have been neces- 
sary to achieve a significant detection. No evidence for variability was found in the Fermi 
integral flux and photon index lightcurves. That there is no evidence for variability from 
PG 1553+113 is consistent with observations of other BL Lacs with Fermi and the VHE in- 



struments. The recently- released Fermi LBAS ( lAbdo et al.ll2009al ) found evidence that BL 
J_j£LCS 3X6 less variable than the other blazar subclass, the flat spectrum radio quasars (FS- 
RQs). Using a simple \ 2 test, 70% of the LBAS FSRQs were found to be variable compared 
to 29% of the LBAS BL Lacs. Also, the fact that approximately 70% of the EGRET-detected 
blazars are not in the LBAS, with comparable flux thresholds, is a further indication that 
high activity in the gamma-ray range is not frequent for a given source. As the number 
of BL Lacs detected in the VHE regime increases, it seems that these objects are not all 
as variable in t he VHE regime as initial observations of Markarian421 and Markarian 501 



suggested (e.g.. iBuckley et al.lll996l iGaidos et al.lll996l . lAharonian et al 



19971. iQuinn et al. 



(I Wakely &: Horan 



1999h . Of the 21 BL Lacs now confirmed to be VHE gamma-ray emitter; _ 
20081 ). 10 (48%) have, to date, shown no evidence for strong variability, 4 (19%) have shown 
marginal evidence for variability (i.e., different flux levels on yearly timescales), and 7 (33%) 
show strong evidence for variability. It should be noted however, that selection effects can 
play a role here: the low duty cycles and small fields of views of VHE instruments limits their 
ability to perform routine monitoring of a large sample of sources and, in many cases, BL Lac 



8 http://tevcat. uchicago.edu; The "Default Catalog", which contains only those sources that have been 
confirmed as VHE emitters, was used (see TeVCat website for more details). 
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Table 1. The flux and spectral parameters from X-ray observations (2-10 keV, unless 
otherwise noted) of PG 1553+113. The spectral parameters for a log-parabola fit are shown 
with the differential X-ray Flux, Fx(E) oc E - " - ^ 10 ^^ (cm -2 s _1 keV -1 ). In cases where no 
value is given for b, the best-fit power-law spectrum is shown with Fx(E) oc E~ a 



'cm 2 s 1 keV 1 ). 



'he abbrevi ated references corre spond to t he following: [DonOl] 



Donato et all fl200lh: [D on05] - iDonato et all fl2005h: [Per05] -IPerlman et all f|2 005): 



[Qst06] - lOsterman et al' l d2Q06h : [Tra07l - iTramacere et all J2007h : [Rei08] - iReimer et al. 



fl2008h : [Mas08] - iMassaro et all hoO$ ). 







Flux (2-10 keV) 




Spectral 




Observation 


(x 10~ n ergcm~ 2 s _1 ) 




Parameters 


Observatory 


Date(s) 


[Ref.] 


a 


b [Ref.] 


ASCA 


1995-08-16 


2.9 [DonOl] 


2.47 


[DonOl] 


BeppoSAX 


1998-02-05 


1.4 [Don05] 


2.17 


0.63 [Mas08] 


XMM-Newton 


2001-09-06 


3.5 [Per05] 


2.09 


[Mas08] 


RXTE 


2003-04-22 - 2003-05-12 a 


0.7 [Ost06] 


2.37 


[Ost06] 




2003-04-22 - 2003-05-28 b 


0.5 [Ost06] 


2.60 


[Ost06] 




2003-04-26 c 


0.3 [Ost06] 


3.19 


[Ost06] 




2003-05-ll d 


1.2 [Ost06] 


2.26 


[Ost06] 


Swift-XRT 


2005-04-20 


2.1 [Tra07] 


2.21 


0.36 [Mas08] 




2005-10-06 


6.9 [Tra07] 


2.14 


0.24 [Mas08] 




2005-10-08 


6.7 [Tra07] 


2.11 


0.23 [Mas08] 


Suzaku 


2006-07-24 - 2006-07-25 


3.5 [Rei08] 


2.19 


0.26° [Rei08] 



a The mean flux and spectral index during the whole campaign. 
b The mean flux and spectral index prior to the flare. 

c The minimum flux recorded during the campaign and the spectral index at that time. 
d The maximum flux recorded during the campaign and the spectral index at that time. 
e The Suzaku spectrum was measured between 0.3-30 keV. 
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Table 2. The differential ffux measured by the Fermi LAT in each energy bin. 



Energy Range 
(GeV) 



Flux 



cm 



- 1 GeV" 



0.20 


- 0.43 


8.63 ±1.56 


X 


10 


-8 


0.43 


- 0.92 


2.57±0.28 


X 


10 


-8 


0.92 


- 1.96 


6.46 ±0.70 


X 


10 


-9 


1.96 


- 4.18 


2.26 ±0.24 


X 


10 


-9 


4.18 


- 8.94 


5.38 ±0.78 


X 


io- 


-10 


8.94 - 


- 19.13 


1.59±0.28 


X 


io- 


-10 


19.13 


- 40.90 


5.52 ±1.08 


X 


io- 


-11 


40.90 - 


- 187.05 


5.37i^t x : 


io- 


12 



^The uncertainty in the highest en- 
ergy bin is found to be asymmetric. 



Table 3. The parameters of the electron distribution in the SSC model, as described in 
the text, used to fit the PG 1553±113 data. The dates of the observations are listed at the 
top of each column. The colour of the line used to show each of the models in Figure H] is 
listed in parenthesis beneath each column's header. For each dataset, the minimal Lorentz 
factor is ^min = 1-00 and we find an equipartition factor of 0.2 (ratio of the energy density 
in the magnetic field to that in the particles of the jet). 



Model 


RXTE 


Swift-XRT 


Suzaku 


Swift-XRT 


No X-rays 


Parameter 


(Blue) 


(Yellow) 


(Green) 


(Magenta) 


(Black) 




Apr/May 2003 


Oct 2005 


Jul 2006 


Mar 2009 




Pi 


1.70 


1.70 


1.70 


1.70 


1.70 


P2 


3.00 


3.00 


2.70 


3.00 


3.00 


P3 


4.10 


4.10 


3.90 


4.10 




Imax x 10 6 


3.16 


3.16 


3.16 


3.16 


0.20 


7i x 10 4 


6.59 


6.59 


6.59 


5.07 


5.71 


72 x 10 4 


7.61 


22.8 


7.61 


6.59 




D tot x 10 54 


4.00 


3.70 


3.93 


4.28 


4.32 
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observations are triggered by high flux states. Fermi, on the other hand, with its large field 
of view and high duty cycle, can sample the lightcurves of sources in a more even fashion, 
in particular over day-scale timebins. It is, however, less sensitive than VHE instruments to 
shorter-term variability of sources with spectral and flux characteristics typical of BL Lacs 
detected to date. Threshold effects, also, can play a role in the detection of variability from 
a source. When the falling edge of the source's SED intersects with the low energy range 
of the instrument sensitivity (this is the case for the FSRQs with Fermi and the BL Lacs 
with the VHE instruments), a slight spectral shift from the source can mimic the effects of 
a large level of variability. 



4.1. Modeling the Intrinsic Emission of PG 1553+113 

In accordance with the weak observed variability, all of the available PG 1553+113 
gamma-ray data were combined in turn with each of the X-ray datasets and, assuming 
a redshift of z=0. 75 as found above, each dataset was modeled with a single-zone syn- 
chrotron self-Compton (SSC) model, that is, a scenario in which one population of electrons 
is responsible for the broadband emission, producing synchrotron radiation in the radio 
to X-ray regime and upscattering these synchrotron photons to produce the gamma-ray 



emission (IBand &: Grindlavl Il985l ; the model employed here is similar to the one used in 
Aharonian et al.l (120091 ) ) . Figure H] shows the result of this fitting. The electrons are pa- 
rameterized as a three-component power-law, dn/d^ cx r y~ Pl (i = 1 — 3), with minimal and 
maximal Lorentz factors ^ m i n and 7 maa; , break Lorentz factors of 71 and 72 and total electron 
number of D tot . For each of the X-ray datasets, the SED was modeled by finding the pa- 
rameters of the electron distribution that provided a good fit to the shape of the low-energy 
component (radio to X-ray), while keeping the same values for the emitting region radius, R, 
the magnetic field strength, B, and the bulk Doppler factor, Sbuik throughout. This fitting 
procedure resulted in values of R = 1.4 x 10 18 cm, B = 0.01 Gauss and Sbuik = 32. Table [3] 
lists the best fit values of the electron parameters for each of the X-ray datasets. A zoom 
of the high-energy component of the SED for the Fermi data with the dataset that is most 
simultaneous (that corresponding to the KVA, Suzaku and VHE data from July 2006), with- 
out the effects of EBL absorption, is shown in Figure [5j It can be seen that, even with the 
presence of the intrinsic curvature that is inherent to the SSC model, most of the curvature 
in the high-energy portion of the spectrum can be accounted for by EBL absorption. 

We found that this simple, one-zone SSC model provided a reasonably good fit to 
the PG 1553+113 SED. By altering only the distribution of the electrons that produce the 
synchrotron emission, a good fit to the overall SED was found for each of the X-ray flux 
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states and, for all of these model realisations the VHE component of the SED did not 
change significantly: the magnitude of the changes in the SED above ~200 GeV are on the 
order of the VHE statistical measurement uncertainties. This consistency of spectral shape 
implies that the gamma-ray flux could remain consistent with the state seen by H.E.S.S. 
and MAGIC during the observations of 2005-2006, even in the presence of the large changes 
in the X-ray flux level that have been detec ted. Such behay i our w as also observed with 



simultaneous datasets from PKS 2155-304 by lAharonian et al.l ( 120091 ) . In a quiescent state 



during these observations, no correlation was found between its X-ray and VHE gamma-ray 
fluxes. During previous observations however, when it was in a flaring state in the gamma- 



ray r egime, there was strong correlation between the X-ray and VHE fluxes (jCostamante 



20081 ). This behaviour is an indication that the hard X-ray flux of BL Lacs can change 
significantly without resulting in detectable activity in the gamma-ray regime, except for at 
the peak of the SED at these energies, i.e., that X-ray variability can be accompanied by 
VHE gamma-ray quiescence and a measurable shift in the spectrum at the peak of the Fermi 
energy regime. In such a scenario, the electrons producing the variable X-ray emission are 
at higher energies than those upscattering the bulk of the synchrotron photons to the VHE 
gamma-ray regime; the scatterings of the variable hardest X-rays are suppressed mostly due 
to Klein-Nishina effects but also because of the decreasing target photon density at these 
energies. This effect is demonstrated for the extreme case in Figure H] where the black dashed 
curve shows the broadband spectrum when the high-energy electron component is omitted. 
It can be seen that the difference between it and the curves in which the highest energy 
electrons are included is still on the order of the statistical measurement uncertainties in the 
VHE regime. 

In the framework of this model, the X-ray and VHE gamma-ray fluxes would be corre- 
lated during gamma-ray flaring states while during the more common gamma-ray quiescent 
states changes in the X-ray flux would not lead to detectable changes in the gamma-ray flux, 
except for at the high-energy peak in the SED. The simple, one-zone SSC model employed 
here allows for such a scenario — accounting for the X-ray flux variations observed histori- 
cally while not requiring detectable changes in the VHE gamma-ray flux state. It provides 
a good fit to the VHE gamma-ray data with each of the low-energy datasets, only one of 
which, drawn in green in Figure HI is quasi-simultaneous (the KVA, Suzaku and all but 7 
hours of the VHE data are from July 2006). It remains to be seen whether such behaviour 
is observed in general from BL Lacs, an investigation that can be undertaken jointly by 
Fermi and the X-ray and VHE observatories. Correlations between the arrival of >100 GeV 
photons in the Fermi data and increasing of the hard X-ray flux, in the absence of VHE 
variability, would be indicative of the proposed scenario. Such effects would be difficult to 
pin down however, due to the low rate of detection of of >100 GeV photons with Fermi. 
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Vercellone et all (120041 ) studied the duty-cycle of gamma-ray blazars by comparing the 



sources detected by EGRET with a sample of radio-selected candidate gamma-ray blazars. 
They found that most blazars have a duty cycle of less than 10%, meaning that they spend 
more than 90% of their time in a non-flaring gamma-ray state. The gamma-ray data pre- 
sented here suggest that PG 1553+113 was in such a state in the Fermi energy range dur- 
ing all of these observations. Given the increase in sensitivity and duty cycle afforded us 
by Fermi, we should expect to detect many more blazars in their quiescent states in the 
>200 MeV energy range. 



5. Conclusion 



We have combined the Fermi data from PG 1553+113 with data from radio through 
VHE gamma rays to study its broadband emission. We demonstrated that a simple, one- 
zone SSC model provides a reasonably good fit to the observational results: It accounts for 
the different X-ray flux states observed while also allowing the gamma-ray data to remain 
approximately constant. More detailed theoretical modeling of the SED, which naturally 
might fit the data points better, is beyond the scope of this paper. 

We have used gamma-ray data spanning four ord e rs of magnitude to seek the EBL 
column density (as parameterized by iFranceschini et all (120081 )) that best fits the measured 
spectrum of PG 1553+113. We find that an EBL integrated to a redshift of z = 0.751q!os 
provides the best fit. Such a high value for the redshift would make PG 1553+113 the 
most distant source to be detected in the VHE regime, an attribute consistent with it being 
the VHE source with the largest spectral break in the gamma-ray regime (T VHE - T HE 



~ 2.7). Three assumptions were made in our redshift calcu 
the EBL behaves in a manner consistent with the model of 



ation the first one being that 



Franceschini et all (120081 ) . Of 
the many EBL mode ls available in the l it eratu re (see Finke &: Razzaque ( 20091 ) for a recent 



comparison), that of IFranceschini et all (120081 ) predicts the lowest level of EBL across the 
energy range of interest here; it provides the minimum level of EBL photons that could exist 
based on known sources alone and therefore, in using it, the redshift derived for PG 1553+113 
can be considered an upper limit. Had one of the other EBL models been used, the best-fit 
redshift would have been lower. Indeed, a more in-depth study of the absorption effects of 
the EBL on the spectrum of PG 1553+133 should consider all of the EBL models available 
in the literature. In this way, the systematic effects of using different realisations of the 
EBL could be determined; the data points provided in Table [2] allow for such a study to be 
undertaken. 



Two additional assumptions were made in our redshift calculation, firstly, that the emis- 
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sion state of PG 1553+113 did not change significantly during the gamma-ray observations, 
which can be interpreted as an indication that it was in a low-flux state state at these en- 
ergies throughout those observa tions. Indeed, the quie scent state is the state in which we 
are most likely to find a blazar (jVercellone et al. 2004 ). Secondly, it was assumed that the 
power-law spectrum measured by Fermi does not suffer from significant intrinsic absorption 
up to 1 TeV. Should there be intrinsic absorption of VHE gamma rays at the source, a lower 
level of EBL absorption would be required to best fit the measured VHE spectrum. We note 
that, since the Fermi data show no evidence for spectral curvature, if it is significant intrinsic 
absorption that is mostly responsible for the sharp break in the gamma-ray spectrum, its 
effects would have to kick in at energies exactly above those accessible to Fermi and, at a 
significant level, in order to produce the sharp break observed in the data. 

This is the first time that the measurement of the complete HE to VHE gamma-ray 
spectrum of a source ha s been used to constrain its r edshift. The value derived here is close 
to the limits de rived by lAharonian et all (120081 ) and lAlbert et all (120071 ) and is higher than 
that derived by iMazin fc Goebell (120071) . Differe nt pa rameterizations of the EBL were used 
in these analyses (those of iKneiske et all (120041 ) and iPrimack et all (120011 )) so this is one 
of the factors contributing to the different results. Additionally, for these previous redshift 
constraints using gamma-ray data, it was necessary to make assumptions about the hardness 
of the intrinsic spectrum in the HE regime. The Fermi measurement of the HE spectrum 
from PG 1553+113 allows us to reduce the number of assumptions used in the constraining 
the redshift. Combining the Fermi measurement with the measured VHE spectra affords us 
complete coverage of the PG 1553+113 SED from 200 MeV to 1.1 TeV. 

Its potentially large distance together with its consistently detectable gamma-ray flux 



fects first proposed by 


Protherod ( 


1986) and, subsequently explored by many authors (e.g., 


Prot 


leroe & Stanev 1993; 


Aharonian et al.lll994: Biller 1995: D'Avezac et al. 


2007: Elviv et al. 


2009 


). If present, pile-up from such cascades is predicted to occur below 100 GeV. Some au- 



thors predict that cascades produce a characteristic spectral bump, which could be detectable 
in the Fermi energy regime if the strength of the extrag alactic magnetic field i n the direc- 
tion o f the source being studied is less than B = 10~ 6 nG (ID'Avezac et alll2007l ). lElyiv et al 



(120091 ) predict that, for magnetic field strengths of £><10" 7 nG, extended emission due to 
the cascading of the source photons should be detectable around extragalactic gamma-ray 
sources by Fermi. 

The possibility that PG 1553+113 could be a particularly distant TeV source provides 
further impetus for IR/optical/UV astronomers to revisit the redshift measurement. A 
direct measurement would be very welcome, and would ultimately settle the question of the 
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accuracy of estimates based on the EBL. The Fermi LAT is continuing to accumulate data 
on PG 1553+113. Subsequent studies of these data will, therefore, allow us to measure its 
spectral shape with greater sensitivity at the highest energies accessible to Fermi 
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Fig. 1. — The differential spectrum from PG 1553+113 as measured by Fermi . The solid 
line shows the fit of a power-law to the overall spectrum derived from all of the data with 
energy E > 200 MeV. The data-points (crosses) indicate the flux measured in each of the 
eight energy bins indicated by the extent of their horizontal lines, when the data in these 
energy ranges were analyzed with the photon index fixed at the value derived from the entire 
dataset. The grey shaded area shows the extent of the Fermi 68% confidence band. 
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Fig. 2. — The lightcurve for PG 1553+113 measured by Fermi between MJD 54683 and 
54883 (2008-08-05 - 2009-02-21). The data are binned in 2-day bins. The top panel shows 
the integral flux, I 2 Da y (E >200MeV), for each 2-day bin. The bottom panel shows the 
power-law photon index, T 2 Day, for each 2-day bin. The integral flux and photon index, with 
their uncertainties, calculated by gtlike for the entire dataset are shown by the shaded 
horizontal bands. For twelve of the timebins, indicated by the triangles in the lower panel, 
the analysis failed to converge with reasonable error bars and therefore, these data were 
excluded from the variability analysis and the points are shown as "x" symbols. 
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Fig. 3. — The spectral energy distribution for the gamma-ray data. The individual datasets 
are described in the text. The Fermi datapoints are shown as black crosses. The EGRET 
upper limit for emission above 100 MeV is shown as a black triangle. The H.E.S.S. data 
combined from May & August 2005 are shown as blue dots and those combined from April 
& July 2006 as blue open circles. The MAGIC data combined from April & May 2005 and 
from January to April 2006 are shown as cyan x's while those from July 2006 are shown as 
cyan open squares. The green solid line shows the power-law fit to the Fermi data, extended 
to higher energies with the level of EBL that best fitted the VHE data, which corresponds 
to the EBL integrated to a redshift of z = 0.75. The green dotted line shows the extension 
of the Fermi power-law to higher energies without absorption. The upper and lower 68% 
uncertainty-contours for the Fermi data are shown as red dashed and blue dash-dotted lines, 
respectively. Each of them were also extended with the level of EBL that best fitted the 
VHE data, which corresponded to a redshift of z — 0.79 for the upper contour and to a 
redshift of z = 0.70 for the lower contour. Their unabsorbed extensions to higher energies 
are shown as red and blue dotted lines. 
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Fig. 4. — The spectral energy distribution for PG 1553+113 fit with an SSC model. The 
fitting procedure is described in the text. The yellow, green, blue and magenta lines are, 
respectively, the SSC model tuned to fit the Swift XRT data, the Suzaku data, the RXTE 
data and the Swift XRT data. The black dashed line is the SSC model when the high- 
est energy electrons are omitted entirely from the fit. The parameters are described in the 
text and are summarised along with the dates of the observations in Table [3j The X-ray 
data have been de-absorbed for a column density of Ng = 3.67x 10~ 2 cm" 2 . Apart from 
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20081 ) observations were made in July 2006 and are thus quasi-simultaneous. They 
are shown as green fil led circles. The Swift-XRT and UVOT data from October 2005 
( iTramacere et al.l 120071 ) are shown as yellow triangles while the magenta open squares are 
Swift-XRT data (Swift observation ID 31368001) from March 2009. The RXTE data from 
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Fig. 5. — A zoom on the high-energy portion of the spectral energy distribution for 
PG 1553+113. The Fermi datapoints are shown by black filled circles; the 2006 H.E.S.S. 
and MAGIC data are shown as green solid squares while the H.E.S.S. and MAGIC data 
from 2005 are shown as grey solid circles. The EGRET upper limit is shown as a black tri- 
angle. The best SSC model fit for the Suzaku X-ray dataset, with EBL absorption applied, 
is shown as a green dashed line. The green dash-dotted line shows the SSC model before 
absorption for the EBL. The bl ack solid line shows the F ermi power-law spectrum absorbed 
for z = 0.75 using the model of iFranceschini et al.l (120081 ) while the dotted black line shows 
the unabsorbed Fermi spectrum. The shaded area shows the Fermi butterfly and the grey 
dashed lines, its unabsorbed extension to higher energies. The grey solid lines show the 
Fermi butterfly absorbed for the best fit redshift to each edge, as discussed in Section [31 
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